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Abstract 

For a class of solutions of the fundamental difference equation of isotropic loop quantum cos- 
mology, the difference equation can be replaced by a differential equation valid for all values of 
the triad variable. The differential equation admits a 'unique' non-singular continuation through 
vanishing triad. A WKB approximation for the solutions leads to an effective continuum Hamil- 
tonian. The effective dynamics is also non-singular (no big bang singularity) and approximates 
the classical dynamics for large volumes. The effective evolution is thus a more reliable model for 
further phenomenological implications of the small volume effects. 
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I. INTRODUCTION 



The singularities of classical general relativity, when specialized to homogeneous, isotropic 
models, manifest as reaching zero physical volume at finite synchronous time in the past. 
This in turn imply unbounded growth of space-time curvature and of matter densities etc 
and signals break down of the evolution equations at finite time in the past. It is widely 
believed that this feature of the classical theory will be modified in a quantum theory of 



gravity and recent development of Loop Quantum Cosmology (LQC) 
corroborate this expectation 0, 1^ . 

The mechanism of 'singularity avoidance' involves replacement of the classical evolu- 
tion equation by a quantum one which is a difference equation [4] thanks to the necessity of 
using holonomy operators in the quantization of the Hamiltonian constraint in LQC. This 
equation exhibits the property that the quantum wave function can be evolved through 
zero volume unambiguously. In addition, the discreteness of the triad operator (having zero 
eigenvalue) necessitates defining inverse triad operator (or inverse scale factor operator) 
indirectly. Thanks to the loop representation on the (non-separable) kinematical Hilbert 
space of LQC 0], these operators get so defined as to have a bounded spectrum implying 
only a bounded growth of curvatures/matter densities. This is true for all allowed values of 



the ambiguity parameters 
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lul. 

The non-separable structure of the Kinematical Hilbert space of LQC however, also im- 
plies a huge set of solutions of the Hamiltonian constraint (a continuous infinity in the 
gravitational sector alone). Presumably, a suitable choice of physical inner product can be 
made to cut down the size of the admissible solutions of the Hamiltonian constraint. A 
choice of inner product however is not yet available. The exploration of Dirac observables 



is also at a preliminary stage 



12|. The general issue of whether or not the non-separable 
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14| . In the present 



kinematical Hilbert space is mandatory, is currently an open issue 
work however we assume the current framework of LQC 

Despite the open issues, it is possible to develop a WKB type semi-classical approximation 
from which an effective continuum Hamiltonian constraint can be deduced This 
at once gives access to the usual classical Hamiltonian methods to construct and analyze 



the underlying difference equation to the Wheeler-DeWitt differential equation followed by 
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the WKB ansatz for its solutions. For large volume corresponding to classical regime, the 
continuum approximation is always available, in fact as a requirement on quantization of the 
Hamiltonian operator. In this regime, the WKB ansatz naturally reproduces the classical 
Hamiltonian as the leading o{hP) term. We would like to extend this method also to small 
volumes. 

The large freedom offered by the non-separable structure of the kinematical Hilbert space 
can be exploited to propose a restriction to those solutions of the Hamiltonian constraint 
for which a continuum approximation is valid for all volumes. One can then develop the 
effective classical Hamiltonian constraint for all volumes and explore its consequences. 

In this work we develop such a picture and in comparison with the usual FRW equa- 
tions identify the effective density and pressure which includes the contributions of quan- 
tum fluctuations of the geometry. Some elementary consequences are also noted. Further 



19|. 



implications for phenomenology are discussed in separate papers 

In section II, we detail the effective Hamiltonian constraint for isotropic models. In 
section HI, we discuss the qualitative features of the corresponding dynamics namely, the 
possibility of 'bounce' solutions as well as solutions that could attempt to 'pass through' 
the zero volume and connect to the oppositely oriented isotropic universe. We discuss what 
features of the quantum evolution are captured by the effective classical evolution. In section 
IV, we summarize our conclusions and outlook. 



II. EFFECTIVE HAMILTONIAN CONSTRAINT 

a 



The Kinematical Hilbert space 
of the densitized triad operator. 



is conveniently described in terms of the eigenstates 



The action of the volume operator on the triad basis states are given by 



v\^^) 



1/^) := 1/^) • 



(1) 



(2) 



In the isotropic context we have two classes to consider, namely spatially fiat and close 
models. The quantization of the corresponding Hamiltonian operators is given in |5|. By 
introducing a parameter ri we can deal both classes together. The values ?7 = and = 1 
will give the fiat and the close models respectively. 



The action of the gravitational Hamiltonian on the triad basis states is then given by 

(e-Mo.|^ + 4^^) - {2 + A^l^^i^)\^) + e^^«^|/x-4/io)) • (3) 

Here, /io is a quantization ambiguity parameter which enters through the fiducial length 
of the loops used in defining the holonomies. It is a real number of the order of 1. Notice 
that the Hamiltonian connects states differing in their labels by ±4yUo- This is a direct 
consequence of the use of holonomy operators which have to be used in the quantization of 
the Hamiltonian operator and is responsible for leading to a difference equation below. 

A general kinematical state |s), in the triad basis has the form 

|s) = \^) ( sum over countable subsets) . (4) 

The Hamiltonian constraint of the classical theory is promoted as a condition to define 
physical states, i.e., 

(^S+^StL)N) = 0. (5) 

In terms of 5^ := e'^^s^ and The Hamiltonian constraint (0) translates into a difference 
equation, 

= ^^l+i^n)S^+'^^lQ — (2 + 4/iQ7^?7)A^s^ + ^^^-4/^0 ■5/^-4/^0 

+ 8ft:7wQ7^p) 'hMs^ , V/i gM (6) 

Af, := |/i + /io|^ - |/^-/^0|^ , ^matlrl;") — H^lf^) ■ 

Hjnifi) is a symbolic eigenvalue and we have assumed that the matter couples to the gravity 
via the metric component and not through the curvature component. In particular, Hm{fJ' = 
0) = 0. There are a few points about the above equation worth noting explicitly. 

Although /i takes all possible real values, the equation connects the coefficients only 
in steps of 4fiQ making it a difference equation for the coefficients. By putting fi := u + 
(4/io)n, n G Z, z/ G [0,4/io), one can see that one has a continuous infinity of independent 
solutions of the difference equation, labelled by u, S'^ := Sj,+4^g„. For each v an infinity 
of coefficients, 5*^, are determined by 2 'initial conditions' since the order of the difference 
equation in terms of these coefficients is 2. Coefficients belonging to different v are mutually 
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decoupled. Since the coefficients and the symbohc eigenvalues Hmifi), both vanish for 
fi = 0, the coefficient So decouples from all other coefficients. 

For large values of /i ^ 4/io {n ^ 1), which correspond to large volume, the coefficients 
become almost constant (up to a common factor of ^/n) and the matter contribution is also 
expected similarly to be almost constant. One then expects the coefficients to vary slowly 
as n is varied. This suggests interpolating these slowly varying sequences of coefficients by 
slowly varying functions of the continuous variable p{n) := ^jipn [l5|. The difference equa- 
tion satisfied by the coefficients then implies a differential equation for the interpolating 
functions which turns out to be independent of 7 and matches with the usual Wheeler- 
DeWitt equation of quantum cosmology. This is referred to as a continuum approximation 



17| . This is of course what one expects if LQC dynamics is to exhibit a semi-classical be- 
havior. While admissibility of continuum approximation is well motivated for large volume, 
one also expects it to be a poor approximation for smaller Planck scale volumes. 

This logic is valid when applied to any one of the solutions 5*^. Thanks to the non- 
separable structure of the Hilbert space, we have an infinity of solutions of the Hamiltonian 
constraint. Although S'^ are uncorrelated for different z/, nothing prevents us from choosing 
them to be suitably correlated. In effect this amounts to viewing s^ themselves as functions 
of the continuous variable fi and stipulating some properties for them. In the absence of a 
physical inner product, we don't have any criteria to select the class of solution. It is then 
useful to study properties of classes of solutions of the Hamiltonian constraint. 

The class that we will concentrate on is the class of slowly varying functions. For 
these we will be able to have a continuum approximation leading to a differential equa- 
tion. Making a WKB approximation for this differential equation, we will read-off the 
effective classical Hamiltonian constraint. In anticipation of making contact with a clas- 
sical description, we will use the dimensionful variable := |7^p/i as the continuous 
variable. Correspondingly, we define po := |7^p/io which provides a convenient scale to 
demarcate different regimes in p. We also use the notation: ilj{p{fi)) := s^. Now the def- 
inition of a slowly varying function is simple: ip{p) is locally slowly varying around q if 
ipiq + Sq) ^ ip{q) + Sq^ + \5q^^^ + ■ ■ ■ with successive terms smaller than the preceding 
terms, for 5q < 4po- It is slowly varying if it is locally slowly varying around every g G M 



17| . Note that even an exponentially rising function can be locally slowly varying if the 
exponent is sufficiently small. 
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To explore the possibility of slowly varying solutions of the difference equation ©, 
consider the difference equation more explicitly. For v G (0, 4/^0)5 putting Sniy) '■= 
Su+ifion, ^ni^) ■= ^u+i^ion and momentarily ignoring the matter term for notational sim- 
plicity, the difference equation (jUI) can be written as, 



(2 + 4/1272^) 



An+li 



Sn{v) (7) 



Its general solution can be written as ^^(z/) = 5'o(z^)p„(z/) + 5'i(z/)cr„(z/), where the p„,o"„ 
are fixed functions of v determined by the same difference equation with the 'initial' 
conditions: po^u) = 1, pi(z^) = and (Jo^u) = 0, cri(z/) = 1 and Sq, Si are arbitrary functions 
of z/ G (0,4/io). (The linearity of the equation means that only the ratio A(z/) := 5'i(z/)/5'o(z/) 
(say) parameterizes the general solution.) 

It is clear that the arbitrary functions allow us to control the variation of within an 
interval of width 4/io. At the integral values of /i/(4/io) corresponding to u = 0, there is a 
consistency condition coming from vanishing of the highest (lowest) order coefficient which 
fixes the ratio of 6*0(0) and 5*1(0). The values of 5^=4^on are fixed (up to overall scaling). 
The slowly varying class of functions will be assumed to approximate these exact values. 
The continuum approximation developed below may not be a good approximation at a finite 
subset of these values corresponding to smaller n. 

With these remarks, we now proceed to derive consequences from the assumption of 

(every where) slowly varying, approximate solutions of the difference equation ©. 

3 

Defining A{p) := (g7£p)2A^ and substituting in terms of slowly varying ip{p) in the 
difference equation (0), leads to the differential equation, 

= Bo{p,po)'ipip) +^poB^{p,po)i!'{p) + 8plB+{p,po)i'"ip) where, (8) 
B±ip,Po) ■■= A{p + Apo) ± A{p - Apo) , and 
Bo{p,po) := A{p + 4po) - (^2 + 144^r/^ A{p) + A{p - 4po) + (sSSk^^ HM 

In the above equation, terms involving higher derivatives of ipip) have been neglected as 
being sub-leading in the context of slowly varying solutions. This is not quite the continuum 
approximation referred to earlier since there is 7 dependence hidden inside po appearing 
explicitly in the coefficients of the differential equation. This is also not quite the Wheeler- 
DeWitt equation since this equation is valid over the entire real line (since p can take negative 
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values corresponding to oppositely oriented triad) while the Wheeler-DeWitt equation using 
the scale factor as independent variable is defined only for half real line. 

From the definitions of the coefficients A, B± q, it is obvious that under p —p (change 
of orientation of the triad), A,B^ and the gravitational part of Bq are all odd while B^ is 
even. For notational convenience we restrict to p > while writing the limiting expressions, 
the expressions for negative p can be obtained from the odd/even properties noted above. 
There are two obvious regimes to explore which are conveniently demarcated by the scale 
Pq, namely, p ^ Po and < p Pq. The corresponding limiting forms for the coefficients 
Bq, B± are easily obtained. One gets, 

1 3 

p>Po : Mp^Po) ~ ^PoVp - -^pI P~^ 

49 3 

: B+{p,Pq) ^ GpoVp- ^Po P~^ 
■ B^{p,po) ^ 12pl + o(po p~5) 

: Boip,po) ^ - 12pIp~I - 432^7^ + 288k^HM (9) 

-t-p -t-p 

1 1 o -3 

p < Po : A{p,po) ^ 3pp^ - gP Po ' 

111 
: B+ip,Po) ~ 3(52 - 3^)pp'^ 

: B4p,po) ^ 24(5i-3i) 
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2 3 

: Bo{p,Po) ^ (3p4)(5^ - 3^ - 2) - 432^7^ + 288k^HM (10) 

{-p -Cp 

Notice that for large volume the explicit po dependence cancels out. This equation corre- 
sponds to the usual continuum approximation which has no dependence on 7 and matches 
with the Wheeler-DeWitt equation in a particular factor ordering. For small volume, the 
Po dependence survives, is non-trivial and the coefficient of the first derivative terms is non- 
zero. In view of the even/odd properties of the coefficients, it follows that the first derivative 
of V'(p) must vanish at p = 0. Further more, even for the fiat model without matter, the 
Bq coefficient is non-zero. Had we extrapolated the Wheeler-DeWitt equation from the 
large volume form, we would not have gotten these terms. Thus the quantum differential 
equation (jSJ agrees with the Wheeler-DeWitt for large volume but differs significantly for 
small volume. 

The small volume form of the equation in fact shows that there are two possible behaviors 
namely ip{p) ~ constant or ipip) diverges as an inverse power of p. Neither of the indicial 
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roots depend on the matter Hamiltonian. The latter solution is not slowly varying and 
the former one implies that the wave function has a non-zero value at p = and the wave 
function can obviously be continued to negative p. Thus the differential equation derived 
for slowly varying functions is both consistent at zero volume and mimics main features of 
the difference equations namely passing through zero volume and matching with Wheeler- 
DeWitt for large volume. 

In the ea.Ue. ,uant.a«o„ of „ models based on point Ko.ono.ies taking values 
in f/(l) representations (separable Hilbert space), the decoupling ofs^ coefficient also implied 
a consistency condition which helped select a unique solution [20] from solutions of the 
Wheeler-DeWitt equation valid at large volume. With the non-separable Hilbert space, 
such a condition can only result from family of coefficients. Nevertheless, one has gotten 
a unique solution (up to normalization) thanks to the slowly varying nature of the solutions. 
It is crucial here that for small p ^ 0, the B_ coefficient has a non- vanishing limit which 
forces the first derivative to vanish at p = 0. (If the single derivative term has been dropped, 
both solutions would have been slowly varying near p = 0.) 

In summary, with the restriction to slowly varying solutions, we have a continuum approx- 
imation (differential equation) valid for all values of the triad. Further more the differential 
equation permits a unique solution (for each matter state) passing through p = 0. 

For future reference we also note that the differential equation admits a 'conserved cur- 
rent'. Taking imaginary part of tp* times the differential equation leads to, 

2poB+J' + B_J = , J := - ^{tfj*)' (11) 
Defining J{p) := f^^{p)J{p) such that J' = determines the function /. Explicitly, 

J{p) = constant (^e^ ^^'^^^ {^*^' - ^'(V'*)'} , J' = ; (12) 
— > constant (p) — ip{4'*y} (p ^ Po) 

constant ^p^"^ j - i^H'*)'} {p < Po) ■ 

We will now go ahead with a WKB type solution and infer an effective classical Hamil- 
tonian. 

Let ip{p) = C{p)e^'^^^\ Substitution of this ansatz in (jSI) leads to a complex differential 
equation involving C{p), $(p). The real and imaginary parts lead to, 

= Bo{p,po) + 4poB4p,po) {{inC)'} + SplB^p^po) {-^ + i^nC)' ^ + (£nC)"|(13) 



= 4po5_(p,po) + 8plB+{p,po) {$" + 2^'{inCy} (14) 

The WKB approximation consists in assuming that the amphtude C is essentially constant 
and the double derivatives of the phase are small compared to the single derivatives. Consider 
the eq. lfT^ under the assumption of almost constant amplitude C{p). Then this equation 
is a Hamilton- Jacobi partial differential equation. These are generally, partial differential 
equations involving only first derivatives with respect to time and position and they always 



have an associated Hamiltonian mechanics 



$'2 



Bo{p,po)-8plB+{p,Po)^ = (15) 

Noting that the Poisson bracket between the triad variable p and the extrinsic curvature 
variable K is ^, we identify $' := and arrive at the effective Hamiltonian as. 



H^^{p,K,<P,p^) := -- 

K 



4po 2po 



1 



288pi] 



{B+{p,po)~2A{p)} 



+K^{p,(j),p^) (16) 

We have multiplied by certain factors so as to get the matter Hamiltonian term appear 
without any pre-factors as in the classical case. The equation p3|) of course implies H^^ = 
and we will interpret this as the modified Hamiltonian constraint equation. The effective 
Hamiltonian is also odd under p —>■ —p modulo the matter term. 

Note that the and the rj dependent terms are o{ip). For large volume, the terms 
enclosed within the braces are vanishingly small and the effective Hamiltonian is indeed 
classical (the matter Hamiltonian receiving corrections from the inverse triad operator also 
goes to the classical form without any ip dependence). For smaller volumes, the quantum 
modifications are present with explicit dependence on ip,po. The approximation used does 
not quite lead to a 'classical limit' due to explicit appearance of £p. 

For small volumes, the equation ()10|) shows that Bo,B^,A all vanish linearly with p 
while i?_ goes to a positive constant. The real and the imaginary parts of the equation, 
equations IT^ . then imply that $' and C both vanish, which is consistent with ip ~ 
constant, as deduced directly from the differential equation (jSI). 

To interpret the effective Hamiltonian constraint as generating the space-time dynamics, 
let us use the identification |p| = a^. One can then obtain the extrinsic curvature K from 
the Hamilton's equation of motion of p as, 

. ^ ^ ^ KdH^ ^ _ B+{p,po) ^^ 
^ ' dt 3 dK 6po 



The large and small volume expressions for the effective Hamiltonian and the extrinsic 
curvature are, 



^efr 
K 

K 




144 



2 - 52 +32 



4 
p 



(18) 
(19) 

(20) 

(p«Po) (21) 



{P > Po) 

+ Hm 



The large volume expressions are the same as for the classical Hamiltonian as expected. 
The small volume expressions are useful in exploring the behavior of the effective dynamics 
close to the classical singularity. 

It is worth expanding on the identification \p\ = a?. The basic variables of LQC are 
first obtained for a general homogeneous Bianchi class A models with the Maurer-Cartan 
forms normalized in the usual manner. Comparing with the metric ansatz then leads to the 
relations \pi\ = 0203 and cyclic. The basic variables of the isotropic models are obtained 
from the above Bianchi ansatz, by putting pi = p, aj = a \fl leading to the identification 
above. Let us denote this scale factor as OBianchi- The corresponding spatial Ricci scalar is 
^i?(aBianchi) = — • ^^c other hand, the standard FRW metric ansatz is so chosen 



-. These two normalizations match 



that the spatial curvature is given by ^R{ap^) = - 
provided aeianchi = ■ To avoid writing the suffixes, we just note that while comparing 
the large volume Hamiltonian constraint with the standard Friedmann equation, one should 
use the replacement a ^ |. This of course is relevant only for the close model. 



III. QUALITATIVE FEATURES OF THE EFFECTIVE DYNAMICS 

In the previous section we derived the effective Hamiltonian constraint (fTBj) . using a con- 
tinuum approximation keeping terms up to second derivatives and using the WKB approxi- 
mation. If we include higher derivative terms that these would give perturbative corrections 
in the large volume. Our focus is however on the small volume regime and leading correc- 
tions which for the matter sector include non-perturbative corrections. For our purposes the 
truncation to second derivatives suffices. 
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For large volume, we already see that the effective Hamiltonian reduces to the classical 

_3 

one to within terms of the order of p 2. We are interested in checking if the effective 
dynamics is non-singular and precisely in what sense. For this it is sufficient to focus on the 
small volume expressions. We will compare the classical Hamiltonian extrapolated to 
small volume and the effective Hamiltonian ((201) • 

Consider first the classical case. As the scale factor goes to zero, the matter density 
diverges either as for pressure-less matter or as for radiation. Correspondingly, if^ 
either goes to a non-negative constant or diverges as . The Hamiltonian constraint then 
implies that K necessarily diverges. As is well known, in both cases the scale factor vanishes 
at a finite value of synchronous time and this of course is the big bang singularity. This also 
suggests a necessary condition for singular evolution: p = should be reachable in finite 
time. Equivalently, if p = is not reachable in finite time, the evolution is non-singular. 

Momentarily, let us assume that for some reason, the matter Hamiltonian vanishes as the 
scale factor goes to zero, then K must remain finite and p = is indeed on the constraint 
surface. Further more p evaluated at p = is also zero implying that p = is a fixed 
point (rather a fixed 'sub- manifold' of the phase space of gravity and matter). The p = 
trajectories of the dynamics are then not accessible in finite time and the evolution is non- 
singular. Clearly (non-) divergence of matter Hamiltonian dictates (non-) singular evolution. 

For generic (non-singular) trajectories, there are two possibilities now. Either (a) p = 
is approached asymptotically as t ^ —00 or (b) the trajectory exhibits a bounce, K = Q 
at a finite, non-zero p. For example, in the case of scalar matter, with LQC modifications 
included, the former is realized for flat models {rj = 0) while the latter is realized for close 



models (r/ = 1) |22|. 



Consider now the quantum case. The matter Hamiltonian is guaranteed to vanish due 
to the inverse volume operator definition. The arguments for p = being fixed point apply, 
iowever, due to the presence of ip term in eq. pn|l . there is a bounce independent of rj 
Iflj ]. The p = is completely decoupled from another trajectories. This is exactly the same 
feature exhibited by the fundamental difference equation. The exact solution = soSf^^ 
completely decouples from all other solutions. The bounce is then a completely generic 
feature of isotropic LQC. But a bounce also provides a minimum volume for the isotropic 
universe whose value is dependent on matter Hamiltonian. Such a natural, generic scale has 
implications for phenomenology as well jl^. 
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Notice that while interpreting the effective Hamiltonian as a (modified) constraint equa- 
tion, we are keeping the kinematics of space-time (a pseudo-Riemannian manifold) intact. 
The modifications imply modification of the dynamical aspects or equivalently of Einstein 
equations. To see the modifications conveniently, let us write the Hamilton's equations in 
a form similar to the usual Raychoudhuri and Friedmann equations in terms of the FRW 
scale factor. 

Write the effective Hamiltonian ()16p in the form, 

H""^ = --{aK^ + r](3) + -u + where, (22) 

K K 

a:=^ , f3:=— , u := ( (5+ - 2A) (23) 
Putting p := ^ ( a is the FRW scale factor and k = IGttG ) and denoting ^ by ' leads to, 

A+A-'4^+M^-'4^) + '^^ (24) 

a"' a"' 3 a* V 9 / 3 

a 8a 



3- . 4 

a 3a* 



-kI [2- — \ H„,- aHL 
a 



(25) 



Comparing with the usual FRW equations, we identify effective perfect fluid density and 
pressure as, 

_ 32 aHm 32 (w 6r]_ f _ lQap\ 
P - 3 + 3Ka4 + ^a^ 9 > ^^^^ 

^.rr 32 a f/ aa'\ 1 32 a f / aa'\ , 

- J(l-"J!^)p-a0+^^l\ (27) 



9 na^ 'IV a J ' ' ' 16 a ^ 
The large and small volume expressions for the effective density and pressure are, 

3 3 1 

for p Po • a ^ -a , P -a , u ^ —-ipa~^; 

peff ^ +8 ^-3^^ _ A^a-s , (28) 

3k 

P^^ - -la~'{aH'J - |-4a-6 ; (29) 

for p < po : a ^ > P ^ ^Po > ^ ^ST^^^" ^^^^^ ' ^^^^^ 

b ■= ^/5-V3 ; 
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P' 



p 



,eff 



,cfr 




2— —a ^Hjn — 



a ^ {Hm + aH'J + 



6(2 - b) 



192 K 




The effective density and pressure receive contributions from the matter sector and also the 
spatial curvature (?7-dependent terms). Apart from these, the homogeneous and isotropic 
quantum fluctuations of the geometry also contribute an effective density and pressure (£p 
terms). While tiny, these are non-zero even for large volume. Notice also that while for 
large volume the spatial curvature does not contribute to the density and pressure, for 
small volume it does. As a by product, we have also obtained the density and pressure for 
matter, directly in terms of the matter Hamiltonian. This is useful because currently LQC 
modifications to the matter sector are incorporated at the level of the Hamiltonian and not 
at the level of an action. Consequently, usual prescription for construction of the stress 
tensor and reading off the density and pressure is not available. These definitions of course 
automatically satisfy the conservation equation: ap' = — 3(P + p). 

Let us consider the vacuum sector, Hm = 0. The more general case of presence of scalar 
field matter is discussed in |l3|. Even for the fiat model (?7 = 0), the o{£p) terms contributes 
to the effective density and pressure. This term in the effective Hamiltonian, Wgg := i^^^, 
is a 'potential' term and will be referred to as the quantum geometry potential. It is actually 
of order after expressing p,po in terms of the /i, po- It is easy to see that the quantum 
geometry potential is odd under p — > —p (since and A are odd) and for p > it is negative- 
definite. Its plot is shown in the figure This immediately implies that in the absence of 
matter (and cosmological constant), all the three terms in the effective Hamiltonian must be 
individually zero which is not possible for the quantum geometry potential. In other words, 
there is no solution space-time. This is in contrast to the purely classical Hamiltonian which 
does give the Minkowski space-time ^ as a solution for the fiat case (77 = 0). This is also 
apparent from the non-zero value of the effective density which prevents the Minkowski 
solution. 

This feature can also be understood in the following manner. The differential equation 
has a unique solution (which is slowly varying every where). This solution is however purely 

^ The Minkowski space-time here refers to Riemann flat metric regardless of its global topology. In the 
cosmological context, the spatial slice is always compact. 
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FIG. 1: The quantum geometry potential. The triad variable p is in units of po while the potential 
Wqg{p) is in units of 2^4^*0 ^''^ 

real and does not admit a WKB form. Consequently, the universe does not admit any 
'classically allowed region' in the WKB sense and thus also does not exhibit a classical 
behavior. Once matter is included, we have again a unique, real solution of the differential 
equation for every matter state. We can now have complex linear combinations admitting 
possibility of regimes of WKB form and corresponding classical behavior. 
It is straight forward to write down a Lagrangian from the H^^ as, 

K 



1 

K 



Po ^2 , f 4 iro/^^N nAf^^^ (32) 



It would be interesting to see if this Lagrangian can be obtained from a specialization of 
a generally covariant action to homogeneous, isotropic metric. 

IV. CONCLUSIONS AND OUTLOOK 

The results of this paper are based on two essential ingredients: the proposal of a contin- 
uum approximation for all volumes exploiting the non-separable nature of the kinematical 
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Hilbert space and the derivation of the effective Hamiltonian via the WKB route. 

The continuum approximation step leads to a differential equation for a (still) quantum 
wave function, ipip)- This equation matches with the Wheeler-DeWitt equation for large 
volume and has important deviations (the first derivative terms) from it at small volumes. 
These deviations allow continuation of the wave function through zero volume just as the 
fundamental difference equation does. For slowly varying solutions, it picks out the 'bound- 
ary' condition ^(0) = 0. Again this is analogous to unique solution (per matter state) 
picked out by the difference equation obtained from the U{1) point holonomies in the ear- 
lier work (20]. Thus the essential features of the fundamental difference equation namely 
non-singular quantum evolution with semi-classical limit are captured by the continuum 
differential equation. 

The effective dynamics specified by the effective Hamiltonian deduced via the WKB 
approximation also reflects these features. The effective dynamics is non-singular, captures 
the decoupling of the = So5^,o exact solution of the difference equation by making the 
classical p = trajectory decouple and reduces to the usual classical dynamics of general 
relativity for large volumes. Since the essential features of quantum dynamics are now 
captured in classical geometrical terms, the effective dynamics is more reliable than the 
usual one and one can now simply work with the effective dynamics to do phenomenology. 
Already, at the qualitative level, one sees that all non-trivial evolutions necessarily show a 
bounce providing a natural scale for say, density perturbations and their power spectra. 

Since the approach draws on the WKB method, a few remarks on the interpretational 
aspects are in order. 

From the continuum quantum dynamics, with the WKB approximation, one obtains a 
Hamilton- Jacobi equation. As a mathematical result, any Hamilton- Jacobi differential equa- 
tion has an associated Hamiltonian mechanics and corresponding trajectories. In our 
context, we are interpreting these trajectories as possible evolutions of the isotropic universe. 
There is at least an implicit implication (or assumption) that a quantum system executing a 
WKB approximable quantum motion physically exhibits a classical motion governed by the 
Hamiltonian associated with the corresponding Hamilton- Jacobi equation. The justification 
for this comes from noting that for large volume we expect the universe to exhibit classical 
behavior and there it is indeed governed by the associated Hamiltonian. For how small 
volumes can we assume this expectation? This question is naturally related to the domain 
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of validity of the WKB approximation. 

One expects the WKB approximation (slow variation of the phase and almost constancy 
of the amplitude of the wave function) to break down closer to the classically indicated 
singularity at zero volume. Noting that the differential equation is local in p and its solutions 
are also local solutions (valid in open intervals in p), we can begin with a WKB approximable 
solution valid in the larger volume regime and attempt to extrapolate it to smaller and 
smaller volumes. All through these extensions, one will have the effective Hamiltonian with 
its associated trajectories which can access the values of p in these intervals. The effective 
Hamiltonian constraint defines a submanifold of the phase space and all trajectories must 
lie on this. The range of configuration space variables (eg p in our case) allowed by the 
submanifold defines 'classically allowed region'. As is well known from the usual examples in 
non-relativistic quantum mechanics, the WKB approximation breaks down at the 'turning 
points'. These correspond to the boundary of classically allowed region which therefore 
demarcates the domain of validity of WKB approximation. Clearly, when such a boundary 
is reached by a trajectory, it must turn back. This is of course the bounce (a = 0). The 
expectation that WKB breaks down at non-zero volume translates into the expectation of 
a bounce occurring at non-zero scale factor. The bounce can thus be understood as the 
smallest volume (or scale factor) down to which one may use the classical framework with 
some justification but below it one must use the quantum framework. More details of the 
bounce picture as well as its genericness are discussed in jiol . 

Further justification comes from other known examples. For example, solutions of the 
Maxwell equations in the eikonal approximation can be understood in terms of the normals 
to the wave fronts which follow null geodesic. The interpretation that this actually reflects 
rectilinear motion of light, may be justified by noting that the Poynting vector (energy flow) 
is also in the same direction as the normals. Likewise, in the context of usual Schrodinger 
equation of particle mechanics, the conserved probability current also points along the nor- 
mals to the wave fronts giving credence to the interpretation that a quantum state of the 
WKB form realizes motion of a particle (or wave packet) governed by the associated Hamil- 
tonian mechanics. In both these examples, further inputs other than the mathematical 
association between Hamilton- Jacobi differential equation and a Hamiltonian system, seem 
to be needed to understand the physical realization of the Hamiltonian system. 

Interestingly, the equation (jH)) does admit a conserved current (fT^ which indeed is pro- 
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portional to the gradient of the WKB phase. Whether this could be used for guessing 
physical inner product vis a vis a probability interpretation remains to be seen. 
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